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Many-body Hamiltonians obtained from first principles generally include all possible non-local 
interactions. But in dynamical mean field theory the non-local interactions are ignored, and only the 
effects of the local interactions are taken into account. The truncation of the non-local interactions is 
a basis dependent approximation. We propose a criterion to construct an appropriate localized basis 
in which the truncation can be carried out. This involves finding a basis in which a functional given 
by the sum of the squares of the local interactions with appropriate weight factors is maximized 
under unitary transformations of basis. We argue that such a localized basis is suitable for the 
application of dynamical mean field theory for calculating material properties from first principles. 
We propose an algorithm which can be used for constructing the localized basis. We test our criterion 
on a toy model and find it satisfactory. 
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I. INTRODUCTION 

In the last decade and a half dynamical mean field the- 
ory (DMFT) has emerged as an important tool for study- 
ing condensed matter systems with strong correlation. 1 
The principal difficulty in understanding these systems is 
the non-perturbative character of such systems, for which 
the physical properties cannot be understood by expand- 
ing various quantities in powers of the interaction. In this 
respect DMFT is a powerful tool for studying problems of 
interacting electrons on a lattice. It is a non-perturbative 
technique which is able to capture fully the local dynam- 
ical correlations in the system. Single site DMFT, as an 
approximation scheme, is controlled in that the result is 
exact in the limit of large coordination numbersi Re- 
cent extensions of DMFT to clusters seem to be rapidly 
convergent for local observables^ Other cluster schemes 
such as cluster perturbation theory and dynamical clus- 
ter approximation are also being used to study problems 
of strong electron correlation^ Recently it has also been 
recognized that DMFT can be used as a powerful tool 
for the realistic computation of properties of materials 
as in the LDA+DMFT schemei^ Indeed results for a 
large variety of materials ranging from Cerium^ Iron and 
Nickel, 8 Plutonium^ and many other oxides have been 
successfully studied with this method starting from first 
principles. 

A common way to utilize DMFT in first principles cal- 
culations is to first derive a Hamiltonian with a kinetic 
energy part and a general short-range interaction part. 
This Hamiltonian, which will be the starting point of 
this paper in Eq. (JTJ, is subsequently studied by DMFT. 
The long-range part of the Coulomb interaction can be 
taken into account by several means. For example, in 
extended DMFT this is done by coupling the electron 
at the impurity site to a bath of bosons whose spec- 
tral function is determined self-consistently»i2iii This is 
equivalent to treating the electrons in the presence of a 
fluctuating electric field (long range interaction). An- 
other possibility is to follow along the lines of Bohm 



and Pines^ In this method, starting with the charged 
electron gas, one performs canonical transformations to 
screen the electrons. In the resulting Hamiltonian the 
excitations are no longer the bare charged electrons, but 
screened neutral quasi-particles. There are various meth- 
ods to obtain the starting Hamiltonian for DMFT. (1) In 
one of the approaches the kinetic energy term is the Kohn 
Sham Hamiltonian of a density functional theory calcu- 
lation written in a local basis set. The interaction terms, 
which can include on-site (Hubbard ) as well as the short 
range part of the Coulomb interaction, is evaluated us- 
ing constrained LDA^ (2) In an alternative procedure, 
as mentioned above, one could start with the electron 
gas Hamiltonian and the periodic potential, and perform 
the Bohm-Pines canonical transformation^ to reduce the 
range of the Coulomb interactions, and then write the 
transformed Hamiltonian in a local basis set. (3) A third 
approach proposed recently^ uses the GW approach to 
obtain the interaction strength. 

The next step is the study of the resulting Hamilto- 
nian with a short-range interaction using DMFT. This 
involves local approximations, and the notion of local- 
ity depends explicitly on the basis set considered. To 
illustrate the point, if we perform an invertible trans- 
formation of the original basis, we merely re-express the 
original Hamiltonian in a new basis, provided we keep all 
the terms in the Hamiltonian. The full electron Green's 
function is obtained by applying the same transforma- 
tion to the creation and destruction operators. But in 
practice, one performs two approximations that explic- 
itly depend on the basis set. The first one is to neglect 
interactions whose range exceeds the cluster size (trun- 
cation) . The second (local approximation) consist of set- 
ting equal to zero the elements of the self energy which 
exceeds that size. These two approximations explicitly 
depend on the definition of locality which is encoded in 
the basis set. In this paper we address only the first is- 
sue, and argue that truncating non-local interactions is 
appropriate when the wave- functions of the basis are well 
localized. As DMFT techniques are beginning to be ap- 
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plied to Hamiltonians with realistic interactions involving 
non-local termsir* there is need for well-defined criteria 
for choosing optimal bases for computations. The pur- 
pose of this paper is to propose one criterion which can 
be used to construct a localized basis for DMFT compu- 
tations. 

The method of choosing a suitable localized basis 
of wave-functions has been studied earlier in quantum 
chemistry and in band structure theory^ The formu- 
lation of the problem consists of two steps. First, one 
identifies a certain group of transformations of the basis 
states, say for example, unitary transformations. Second, 
one identifies a criterion that picks out one basis out of 
all possible choices that are connected by the transforma- 
tions. The criterion is a basis dependent quantity, and 
therefore is a functional in the space of the transforma- 
tions. It is a measure of the amount of localization of the 
wave-functions in a given basis. For example, in quan- 
tum chemistry "energy localized molecular orbitals" have 
been studied^ These are obtained by maximizing under 
unitary transformations a functional given by the sum of 
the Coulomb self-interaction of the orbitals. Similarly, 
for band structure calculations the use of "maximally- 
localized" Wannier functions has been proposed^ The 
idea is to exploit the freedom that is present in the choice 
of the phases of the Bloch orbitals. With a given set 
of Bloch orbitals one can define a new set by a unitary 
transformation. From each such set of Bloch orbitals one 
can obtain a corresponding set of Wannier functions by 
Fourier transformation. The maximally-localized Wan- 
nier functions are obtained by minimizing the spread 
functional, which is the sum of the second moments of 
the Wannier functions, in the space of unitary transfor- 
mations. More recently, the construction of localized ba- 
sis states has been extended to include non-orthogonal 
molecular orbitals. 17 

The rest of the paper is organized as follows. In section 
lUwe identify a criterion for choosing a basis suitable for 
DMFT. We construct a functional which is maximum in 
the preferred basis. We discuss the properties of such a 
basis by studying linear variations of the functional under 
unitary transformations. We also propose a method for 
constructing the preferred basis. In section lTrfl we test the 
criterion on a Hamiltonian whose interaction is taken to 
be simple but non-trivial. We find that the criterion and 
the associated functional is well-behaved. In conclusion, 
we summarize our main results. 



II. LOCALIZED BASIS FOR DMFT 

To keep the discussion general, in the following we for- 
mulate the problem in a basis which is non-orthogonal. 
For this purpose we consider a system of interacting elec- 
trons on a lattice whose Hamiltonian is expressed in a 
basis of atomic orbitals. The single particle states are 
denoted by 4> a (r — R„) = (r\na), where a is a symmetry 
related index (say, orbital) and R„ is a lattice position. 



We suppose there are m orbitals per site such that the 
index a = 1, • • • , to, and there are N lattice sites with the 
index n = 0, • • • , N — 1. We also impose periodic bound- 
ary condition |n, a) = \n + N, a). The states defining the 
basis, unlike those in a Wannier basis, are not orthogo- 
nal. We denote the overlap between any two states by 
O a f3(n — to) = (na\m[3). The second quantized many- 
body Hamiltonian can be written as 

n = j2 tn a p c iAp+ E ^ ! <o c i/MCi, 7 - (i) 

nm nmlk 

We assume that the matrix elements i"™ = 
(na\Ho\m/3) for the non- interacting part, and V^pty = 
(na, m(3\V\l"/, k5) for the interacting part are known 
from first principles studies such as band structure cal- 
culations. It is useful to bear in mind that the anti- 
commutation relation between the creation and anni- 
hilation operators in a non-orthogonal basis is given 
by {c^q i c m ^} = 0~p(n - m)w& We now consider 
an invertible transformation of the single particle basis 
that preserves the lattice translation invariance, \na) — > 
\n'a') — J^mpTpa'i'm — n')\m(3). Expressed in the new 
basis the Hamiltonian, say H', has the same form as in 
Eq. Q except with all the indices primed. We know that 
TC = Ti., since it is the same operator expressed in two 
different bases. However, when we truncate all the non- 
local interactions, we deal with a model Hamiltonian of 
the form 

<i_y \ ^ i)im f , \ ^ -trnnnn \ f 

a/3 afl-,5 

(2) 

But the process of truncation is a basis dependent step. 
If we perform the truncation in the new basis, i.e., on 
TL' , the resulting new truncated Hamiltonian H' tI ^ Htr- 
This observation implies that ignoring non-local interac- 
tions is a good approximation only if the single particle 
basis is sufficiently localized. In the following we develop 
a systematic criterion for constructing such a basis. 

Here we consider only unitary transformations of ba- 
sis. Later we comment about the possibility of extending 
the scheme to include non-unitary invertible transforma- 
tions as well. We start from an initial basis {Ina)}, and 
consider unitary transformations 

\na) — » \n a) = U\na) = Up a {m — n)\mf3) (3) 

to new basis states {|n'a')}. In order to find a criterion to 
choose the most localized basis among the possible bases 
{|n'a')}, we first identify a quantity which is invariant 
under unitary transformations. The trace of any oper- 
ator has this property. Since we are concerned about 
truncating the interacting part of the Hamiltonian, wc 
consider the trace of the square of the interaction opera- 
tor. In terms of the overlap matrix and the interactions 
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expressed in the {|na)} basis this is given by 

I = Tv(V 2 ) = 0-^n-m)0-^l-k)0-^r-s) 

x 0${p-q)VffV$? a . (4) 

Here, and in the rest of the paper we adopt the conven- 
tion that repeated indices are summed. The invariant 
defined above has two basis dependent parts, namely, 
terms that involve only the local interactions and those 
involving non-local interactions. Keeping only the local 
interactions in a given basis, we define the "local inter- 
action functional". For example, in the basis {|na)} the 
functional has the value 

F[{\na)}} = O^0)O^(0)O;;(0)O-J(0)^F™. 

To elucidate the structure of the functional we first note 
that the overlap matrix remains unchanged under unitary 
transformations, i.e., 

(n'a'\m' P'} = O a ipi(n' — to') = (na\mf3) = O a p(n — to). 

Next, the transformation of the interaction terms is given 
by 



* ocf3 S-y 



V$F'$ = K a (r-n)U* pf3 (s-m)V™% 
x U m (p - l)U vS (q - k). 



(6) 



In terms of the unitary transformations the local inter- 
action functional can be written as 

F[{\nW)}] 

= o-} , (0)0-), (0)0;,;, (0)O jl, (0)V^P v , a , V™%, 
= '0-1(0)0-1(0)0^(0)0^(0) 



x u; f3 (r)u;: 5 (s)v^zu Ta (p)u Xv ( q ) 

x [u: p (n)U; u ( m )V^JfUg a (l)U^(k)] 



(J) 



The inverse of the overlap matrix enters as weight factor, 
and the interaction terms in the starting basis 
serve as parameters of the functional. The desired basis is 
the one in which the functional is maximum in the space 
of unitary transformations. This criterion also implies 
that, in the chosen basis, the part of the invariant I that 
contains non-local interactions is minimized. 

In order to study the property of the preferred basis 
we consider infinitesimal unitary transformation given by 
U = e leH , where H is hermitian and e is a small parame- 
ter. The action of H on the single particle wave- functions 
is given by H\na) = Hp a (m — n)|TO/3), such that 

U a p(n-m) = S a/3 S nm + (ie)H af3 (n - m) 
(it? 



2! 



-H aj (n — l)Hjfj(l — m) + 



The hermiticity of H implies that 

[{na\H\m(3}]* = (m/3\H\na) 



H, 



7(3 (7 - m)O ia (l - n) 



Of3 7 (m — l)H ia (l — n). 



For a lattice of N sites with periodic boundary condition 
and to orbitals per site, we note that the transformation 
matrix H has Nm 2 real independent parameters. In the 
following we assume that V^ri,^) = V r (r2,ri), so that 
^0/357' = Vfiaryt- For the convenience of notation we 
define the quantity 

L^(t) ee O^(0)O^(0)O^(0)O^(0)^V> ^. (8) 

To 0(e) the variation of the functional can be written as 

SF = (-^e)[L^(t)H; a (t) - L ap (t)H, a (t)] 

= (~4ie) [0-£(m - n)L* af3 (N - m)O ap (n - t) 

-L^(t)}H^(t). (9) 

We define 



A afi (t) ee L^(N-t)-0-^(m-n)L* al3 (N-m)O af ,(n+t), 

(10) 

and we note that A is anti-hcrmitian, i.e., 

0/3 7 (to - l)A ja (l -n) = ~A* p (l - m)O ia (l - n). (11) 

The condition for the functional F to have a local maxi- 
mum is 

SF 



6H^(t) 



A all (N-t)=0. 



(12) 



The above anti-hermitian condition has to be satisfied by 
the preferred basis. In other words, the preferred basis is 
the one in which L api (t) is hermitian. The above condi- 
tion gives Nm 2 real independent equations, which is the 
same as the number of real independent parameters in 
the transformation matrix H. 

The following is a simple ansatz for maximizing F by 
successive unitary transformations. We start with an ini- 
tial basis {|na}}, and we calculate A prT (t) using Eqs. © 
and (|10f> . We then change the basis using the transfor- 
mation 



H^(t) = iA^t), 



(13) 



and iterate this procedure until the condition for the 
maximum is achieved. We assert (proved in the ap- 
pendix) that with this ansatz, to 0(e) 



SF 



-(Ae)A a ^(N-t)A^(t)>0. 



(14) 



This ensures that with successive transformations the 
value of the functional increases (provided e is small 
enough) until it reaches a local maximum. 

Our method for constructing a localized basis is simi- 
lar to Wegner's flow equation 1 — and Glazek and Wilson's 
similarity renormalization scheme^ approaches where 
small unitary transformations are used to reduce non- 
diagonal elements of a Hamiltonian. In our case, we use 
unitary transformations to reduce the strength of non- 
local interactions and increase the value of the functional 
F. 
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III. DISCUSSION AND CONCLUSION 

We test our criterion on a lattice Hamiltonian with an 
interaction of the form 



Hint = ^2 Uijkicjc 



CkCl, 



(15) 



ijkl 



where 



U m = £y (i - R)</>*0' - R)0(* - n-W - R). (16) 

R 

Here <j>{i — R) = ^(R, — R) is an orbital on lattice site 
Ri, which defines a one-particle starting basis {4>(n— R)}. 
The interaction described above is local in configuration 
space, but in the basis of {(f>(n — R)} it has non-local 
terms as well. For simplicity we consider one orbital per 
site, and assume that the starting basis is orthonormal. 
In this basis the local interaction functional is given by 



1 



F = jr E I0(«-R)l 4 l^-R') 

S ra.R.R' 



(17) 



where N s is the number of lattice sites. We note that 
the above functional has a form similar to the (square 
of the) inverse participation ratio that is studied in the 
context of Anderson localization^ We perform a unitary 
transformation of the basis of the form <j){n — R) — > (j)(n — 
R) + 8<f>(n — R), where 

5<f>(n - R) = (ie) ^ W(m — R) + £>(e 2 ). 



The coefficients h mn are hermitian such that h* mn — h nm . 
The variation of the functional F can be written as 

AC 

SF = —r \<P{n - R)| 2 {(t>*{n - B)5(j){n - R) + h.c.} , 



n.R 



where C = ^ R |0(n — R)| 4 is a site independent con- 
stant. Using the above form for 5(j> and the unitarity of 
the transformation, we get to 0(e) 



SF 



4ieC 
Tf7 



(m - R)0(ra - R) 



,R 



|0(m-R)| 2 - |<?!>(n-R)| 2 



The extrema of the functional is given by 



(18) 



R 



0*(m-R)0(n-R) |</>(m - R)| - \<j>{n - R)|' 



= 0, 
(19) 

for all sites (n, m). By inspection, there are two solu- 
tions to the above equation. (1) \4>(n — R)| 2 = l/A/" s ,Vn, 
which is the limit of delocalized states, for which F = 
1/A/" 2 (minimum). (2) \<f>(n — R)| 2 = S n r, which is the 
limit of localized states. In this case F = 1 (maximum), 



and the interaction is entirely on-site. Starting with the 
original basis {<j>{n — R)}, and the ansatz 



h nm = i 22 <i>*{n - H)4>(m - R) 



R 



x [j</>(n - R)| 2 - \<p(m - R) 



(20) 



we get to 0(e) 



SF 



4eC 



J2 |0*(n-R)</)(n-R)| 2 



n.m.R, 



(m-R)| 2 -|0(n-R)| 2 



> 0. 



(21) 



Provided e is chosen small enough (to justify the neglect 
of higher order variations), with the above ansatz it is 
possible to increase the value of F with successive uni- 
tary transformations until the limit of localization is at- 
tained. This simple example illustrates how the local 
interaction functional can be used to construct a basis of 
localized one-particle states. For interactions which are 
more complicated and realistic, it is unlikely that uni- 
tary transformations can make the interactions entirely 
on-site. However, the strengths of the non-local terms 
can be reduced (quantitatively defined by maximization 
of functional F) in a more localized basis. 

Two more comments are of relevance. First, our crite- 
rion ignores the non-interacting part of the Hamiltonian. 
If one starts with nearest neighbour hopping in the origi- 
nal basis, in the localized basis the hopping will be more 
complicated. But the point of view adopted here is that 
the non-interacting part can still be solved exactly. Sec- 
ond, in this paper we consider only unitary transforma- 
tions of basis. This implies that one maximizes the local 
interaction functional within a family of bases with the 
same overlap matrix (say, orthonormal bases, if the orig- 
inal basis is orthonormal). In principle one could probe 
for bases with different overlap matrices by general in- 
vertible transformations. Such a group is non-compact 
and one needs to impose constraints such that the func- 
tional is bounded from above. One possible constraint 
can be imposed in terms of the singular value decompo- 
sition of the transformation matrix, say, the ratio of the 
maximum and the minimum singular values be within a 
specified bound. 

In conclusion, we propose a criterion for constructing 
a localized single particle basis where non-local interac- 
tions can be truncated. Such a basis is appropriate for 
using DMFT for the calculation of material properties. 
We suggest a simple algorithm by which the construc- 
tion of the localized basis can be carried out. By testing 
the criterion on a toy Hamiltonian we conclude that the 
criterion and the associated functional is well-behaved. 

We thank D. Vanderbilt, S. Savrasov, V. Oudovcnko, 
and H. Jeschke for stimulating discussions and useful sug- 
gestions. 



5 



APPENDIX A 



Using the above relation and Eq. (|llf> one can show that 



In this appendix we prove the assertion in Eq. I|14l) . 
First, if the basis is orthonormal to begin with, i.e., 
O a p{n — m) = 5 a p8 nm , it is easy to see that 

Apr® = L, a {N - 1) - L*^(t) = -A*^(t). (Al) 

Then, 5F=(Ae)\A all {N - t)\ 2 > 0. 

If the basis {|na)} is non-orthogonal, we assume there 
exists an orthonormal basis {lar}}} (say, a Wannier ba- 
sis) to which it is related by |ar)) = S(n, a; ar)\na) and 
((ar| = (na\S(n, a; a, r)*. One can show that 



SF 



= (4e)0- /s 1 (m - n)A* a p{m)O ail {n - t)A^(t) 

^ k; t, n)A^ ia (t - m)S(m, a; a, t) 

tm 

(A3) 



N 
> 



4e x - 
AT 2^ 



Q l(n — m) = ,5(n, a; a, r)S(m, /3; a, r)* 



(A2) 
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